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HIGHLIGHTS 


►  A  thermal  fluid  electrochemistry  model  and  a  thermal  mechanical  model  are  presented. 

►  Current  collection  method  significantly  affects  the  output  of  mtSOFC. 

►  Effects  of  Ni  content  on  electrochemical  and  mechanical  performance  are  critically  examined. 

►  Area  current  outputs  of  planar  and  micro-tubular  SOFCs  are  found  to  be  comparable. 
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A  comprehensive  thermal  fluid  electrochemistry  model  for  the  analysis  of  the  electrochemical  perfor¬ 
mance  and  a  thermal  mechanical  model  for  the  analysis  of  the  mechanical  behavior  of  micro-tubular 
solid  oxide  fuel  cell  (mtSOFC)  are  presented.  Material  properties  for  the  models  are  determined  by  the 
available  theory  and  experiment  and  linked  to  the  compositions  and  microstructures  of  the  materials. 
Good  agreement  between  the  theoretical  and  experimental  I-V  relations  is  obtained.  The  mode  is  used 
to  examine  the  effects  of  various  parameters  on  the  electrochemical  and  mechanical  performance  of 
mtSOFC.  Collecting  current  from  both  sides  of  the  anode  is  found  to  significantly  increase  the  cell  output. 
Increasing  the  Ni  content  or  reducing  the  Ni  particle  size  is  generally  helpful  for  improving  the  elec¬ 
trochemical  performance,  but  the  increased  Ni  content  reduces  the  mechanical  stability.  The  suitable  Ni 
content  is  thus  obtained.  The  LSM  content  is  inconsequential  to  the  mechanical  stability  and  should  be 
determined  by  achieving  high  electrochemical  performance.  Properly  designed  mtSOFCs  are  shown  to 
provide  high  current  outputs.  Therefore,  mtSOFC  is  a  promising  technology  with  both  the  benefits  of 
planar  SOFC  for  high  current  density  and  tubular  SOFC  for  thermal  cycling  endurance. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  attractive  power  production  de¬ 
vices  due  to  their  high  energy  efficiency  and  fuel  flexibility  as  well  as 
low  pollutant  emission  [1].  Micro-tubular  SOFC  (mtSOFC)  with  tub¬ 
ular  diameter  under  a  few  millimeters  shows  high  performance  on 
the  thermal  shock  resistance,  volumetric  power  density,  fast  startup 
and  thermal  cycling  and  has  attracted  increased  attentions  [2]. 

The  basic  principles  governing  the  choices  of  SOFC  materials 
have  a  lot  in  common  for  the  tubular,  planar  and  micro-tubular 
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designs.  For  example,  the  electrodes  should  have  high  conductiv¬ 
ities  for  the  charge  transport,  high  catalytic  activities  for  the 
required  chemical  and  electrochemical  reactions,  adequate  poros¬ 
ity  for  the  gas  diffusion,  good  chemical  and  mechanical  compati¬ 
bilities  with  the  electrolyte  and  interconnect  materials  for  the  long¬ 
term  operating  stability  [3].  Therefore,  material  compositions  and 
microstructures  suitable  for  tubular  SOFCs  (tSOFCs)  and  planar 
SOFCs  (pSOFCs)  are  often  the  right  choices  for  mtSOFCs.  However, 
the  mtSOFC  design  has  its  own  distinct  characteristics  that  should 
be  examined  to  realize  its  full  potential.  There  is  a  strong  need  for 
in-depth  analysis  of  mtSOFC  due  to  the  fact  that,  as  a  relatively  new 
design  of  SOFC,  mtSOFC  is  less  studied  and  understood  and  there  is 
a  large  room  for  its  improvement.  As  the  experimental  testing  is 
expensive  and  time  consuming,  its  ability  to  study  the  effects  of 
different  material  parameters  and  operating  conditions  on  the 
performance  of  mtSOFC  is  severely  limited.  Mathematical  models 
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incorporating  the  known  physics  and  behaviors  of  mtSOFC  mate¬ 
rials  to  predict  their  performance  are  important  tools  for  the  un¬ 
derstanding  and  technological  improvement  of  mtSOFCs. 

Numerous  mathematical  models  have  been  built  to  study  the 
effect  of  microstructure  on  the  electrochemical  performance  [4-12]. 
De  Caluwe  et  al.  [10]  developed  a  one-dimensional  button-cell 
model  to  examine  the  influence  of  anode  micro  structure  on  the 
electrochemical  performance.  Joan  et  al.  performed  numerical  ana¬ 
lyses  on  the  optimization  of  cathode  micro  structures  for  anode- 
supported  pSOFCs  [11].  Several  thermal  stress  models  with  varying 
degrees  of  sophistication  [13-24]  are  also  proposed  to  study  the 
mechanical  properties  of  SOFCs.  Yakabe  et  al.  [17-19]  evaluated  the 
residual  thermal  stresses  in  the  electrolyte  of  anode-supported 
pSOFCs  by  a  combined  experimental  and  numerical  analysis.  Seli- 
movic  et  al.  [20]  studied  the  thermal  stress  of  pSOFC  in  both  the 
steady  and  transient  states.  Weil  et  al.  [21  ]  performed  a  finite  ele¬ 
ment  analysis  on  the  magnitudes  of  thermally  induced  stress,  strain 
and  part  deflection  in  the  cell,  seal  and  window  frame  components 
under  uniform  heating  and  cooling  conditions.  Nakajo  et  al.  [22,23] 
used  thermo-electrochemical  models  to  study  the  mechanical  issues 
and  the  current-voltage  ( I-V)  characterization  of  pSOFC  stacks 
during  assembly,  heat-up,  dynamic  operation,  load  shutdown  and 
cool-down  phases.  Recently,  Serincan  et  al.  developed  a  multi¬ 
physics  model  to  investigate  the  effects  of  different  factors  on  the 
thermal  stresses  in  a  mtSOFC  with  the  thermal  expansion  co¬ 
efficients  of  materials  treated  as  constants  [25]. 

Though  there  is  no  shortage  of  research  on  the  numerical 
modeling  of  SOFCs,  the  theoretical  effort  focusing  on  the  effects  of 
electrode  composition  and  microstructure  on  both  the  electro¬ 
chemical  performance  and  mechanical  property  of  mtSOFC  is  rel¬ 
atively  rare  in  the  open  literature.  Flowever,  such  research  effort  is 
quite  important  for  the  understanding  and  development  of  mtSOFC. 
For  example,  a  typical  mtSOFC  consists  of  a  composite  anode  sup¬ 
port,  a  thin  electrolyte  layer  and  a  thin  composite  cathode  layer.  The 
composite  electrode  is  a  mixture  of  electron  conducting  material, 
such  as  Ni  for  the  anode  or  La(Sr)Mn03  (LSM)  for  the  cathode,  and 
the  electrolyte-material,  such  as  yttria-stabilized  zirconia  (YSZ)  [3]. 
An  Ni  content  close  to  40  vol.%  in  the  Ni-YSZ  cermet  may  easily 
produce  an  area  power  density  that  is  higher  than  its  tSOFC  coun¬ 
terpart,  but  substantially  lower  than  that  of  the  pSOFC  design  [3,26]. 
Increasing  the  Ni  content  to  about  50%  may  significantly  increase 
the  power  density  [3],  but  the  increased  Ni  content  induces  me¬ 
chanical  flaw  [27].  Therefore,  both  the  electrochemical  and  me¬ 
chanical  properties  and  their  interplay  should  be  taken  into  account 
in  order  to  properly  balance  their  different  needs  and  improve  the 
mtSOFC  technology  for  practical  applications. 

This  work  describes  a  detailed  multi-physics  model  for  the  nu¬ 
merical  simulation  of  mtSOFCs.  The  physical  properties  of  com¬ 
posite  materials  are  determined  by  the  available  theories  and 
experimental  data.  The  model  is  used  to  examine  the  electro¬ 
chemical  and  mechanical  performances  of  mtOFCs  with  the  rep¬ 
resentative  Ni/YSZ/LSM  material  set.  Effects  of  the  material 
microstructures  such  as  Ni  and  LSM  contents  and  their  particle 
sizes  on  the  electrochemical  performance  and  thermal  mechanical 
behavior  of  mtSOFC  are  systematically  examined.  The  desirable 
material  compositions  and  microstructures  are  thus  deduced. 

2.  Theoretical  method 

There  are  two  main  components  in  our  mathematical  model: 
a  thermal  fluid  electrochemistry  model  for  the  analysis  of  the 
electrochemical  performances  and  a  thermal  mechanical  model  for 
the  analysis  of  the  mechanical  behaviors  of  mtSOFCs.  The  micro¬ 
structure  based  multi-physics  electrochemistry  model  and  the 
structural  mechanics  model  are  solved  sequentially.  The  thermal 


stress  distributions  in  mtSOFC  are  evaluated  with  the  temperature 
fields  obtained  from  the  thermal  fluid  electrochemistry  model.  The 
thermal  stress  distributions  and  mechanical  properties  are  used  by 
the  structural  mechanics  model  to  calculate  the  probability  of 
mechanical  failure  of  the  mtSOFC  materials.  The  electrochemistry 
model  and  the  structural  mechanics  model  may  be  assembled  from 
various  literatures  [7,14,28-31].  The  precise  details  of  the  models 
used  in  this  work  are  described  in  the  following. 

2.1.  Electrochemistry  model 


The  multi-physics  electrochemistry  model  involves  coupled 
electrochemistry,  thermal,  fluid  and  electrical  transport  equations 
with  microstructure-based  properties.  The  basic  model  parameters 
for  the  electrochemical,  electrical  and  thermal  properties  are 
determined  by  the  experimental  data  [32].  The  effects  of  the  elec¬ 
trode  composition  and  microstructure  on  the  material  properties 
are  considered  by  the  microstructure  based  effective  property 
theory  [33].  Effects  of  the  electrode  composition  on  the  electro¬ 
chemical  performance  of  mtSOFC  may  then  be  revealed. 

The  overall  performance  of  SOFC  depends  on  the  operating  cell 
voltage  and  the  output  current.  The  cell  operating  voltage,  Vceii.  is 
lower  than  the  open  circuit  potential  (Nernst  potential)  due  to  various 
polarizations  such  as  activation  polarization,  ohmic  polarization  and 
concentration  polarization,  and  may  be  formally  expressed  as  [30]: 

^cell  =  E(J  ~  Pohm  —  Vcon  ~  Pact *  0) 


where  7jact  is  the  activation  overpotential  for  the  electrochemical 
reactions,  pCon  is  the  concentration  polarization  due  to  the  resist¬ 
ance  for  the  gas  transport  in  porous  electrodes,  p0hm  is  the  ohmic 
overpotential  and  E0  is  the  open  circuit  potential.  The  open  circuit 
potential  of  the  fuel  and  oxidant  used  may  be  calculated  as: 


E0 


A  G° 
1 2F~ 


Ph2(po2)°' 

PS 


(2) 


where  AG°  is  the  free  energy  change  of  the  reaction 
H2  +  0.5O2  ^H20  at  the  working  temperature  and  when  the  par¬ 
tial  pressures  of  all  the  reacting  gas  species  are  1  atm  [34]. 

The  concentration  polarization  (pCon)  is  the  sum  of  the  concen¬ 
tration  polarizations  of  anode,  p£on,  and  cathode,  p£on.  They  may  be 
calculated  as: 


Vco 


=  Sn| 

^Ph2Ph2o 

2F  ' 

^Ph2Ph2o, 

/  n  \  0.5 

JSnl 

2F  ' 

, 

l PoJ 

(3) 

(4) 


where  Ph2(Ph2o)  and  Ph2(Ph2o)  denote  the  partial  pressure  of  H2 
(H20)  at  the  anode  three  phase  boundary  (TPB)  and  the  anode-fuel 
channel  interface,  respectively.  p02  and  p^2  are  respectively  the  par¬ 
tial  pressure  of  02  at  the  cathode  TPB  and  the  cathode-air  channel 
interface. 

The  ohmic  overpotential  includes  the  electronic  and  ionic  ohmic 
overpotentials  of  the  anode,  cathode  and  electrolyte  and  the  con¬ 
tact  overpotential  due  to  contact  ohmic  resistance  at  the  interface 
of  different  cell  components,  which  may  be  written  as: 


Pohm  =  ^ohm  +  ^ohm  +  ^ohm  +  ^ohm  +  ^ohm  +  PASR-  (5) 

Here  and  are  the  electronic  ohmic  overpotentials  of 
the  anode  and  the  cathode,  respectively.  Poiim  and  Pohm  are 
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the  ionic  ohmic  overpotentials  of  the  anode,  cathode  and  electro¬ 
lyte,  respectively.  ijasr  is  the  contact  overpotential  and  may  be 
calculated  as  [30]: 

^ASR  =jmR  ASR>  (6) 

where  j  is  the  local  current  density  at  the  interface  and  Pasr  is  the 
area  specific  contact  resistance. 


2.1.1.  Governing  equations  for  charge  transport 

Electronic  and  ionic  current  densities  in  the  electrodes  and 
electrolyte  are  governed  by  the  charge  continuity  equations  and  the 
Ohm’s  law  [7]: 


V-  i  el  =  V-(-<7fv<pel) 


jTPB^TPB.eff  In  cathode 
_jTPB^TPB,eff  In  anode 


f  -iTPB>TPB  eff  In  cathode 

V-  i  io  =  V-(-(Tf0ffVWo)  =  <J  0  Inelectrolyte  (8) 

{  JTPB^TPB.eff  In  anode 


where  i  el(  i  io)  is  the  vector  of  electronic  (ionic)  current  density, 
°eif(°1of)  *s  effective  electronic  (ionic)  conductivity,  and  <pe\  (<pi0) 
is  the  local  electronic  (ionic)  potential.  ATpB,eff  is  the  effective  vol¬ 
ume  specific  TPB  length,  j’tpb  is  the  TPB  length  specific  transfer 
current  density  and  is  defined  as  positive  when  the  current  flows 
from  the  ionic  phase  to  the  electronic  phase,  j'tpb  is  a  nonlinear 
function  of  the  local  activation  polarization,  i;act,  i.e.,  the  voltage 
sacrificed  to  overcome  the  activation  barriers  associated  with  the 
electrochemical  reaction.  The  relationship  between  j’tpb  and  ^ac t 
may  be  described  by  the  Butler- Volmer  equation  [35]: 


Jtpb  -  Jo 


exp 


2  OfF 

~rT 


-  exp 


RT  Vact) 
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where  af  (ft)  is  the  forward  (reverse)  reaction  symmetric  factor.  j0  is 
the  local  exchange  transfer  current  per  unit  TPB  length  which  is 
estimated  as  [36]: 


Jo  =  Jo,refexP 


(10) 


Jo  =  Jo,refexP 


%  ( 1  1  A  ( Ph2Ph2o\ 

r  v  WJUww 


(11) 


where  T  is  the  temperature  and  R  is  the  universal  gas  constant. 
Jo  ref  (Jo  ref)  an°de  (cathode)  exchange  transfer  current  density 
at  the  reference  temperature,  Tref,  and  is  treated  as  an  adjustable 
parameter  for  fitting  the  experimental  result.  EH2(£o2)  *s  the  acti¬ 
vation  energy  for  the  anode  (cathode)  electrochemical  reaction. 

The  anode  and  cathode  activation  polarizations,  rjlct  and  r]lcV 
for  Eq.  (9)  may  be  expressed  as, 

Pact  =  'Pel  -  <Pio  —  P?on-  (12) 


Pact  —  <Pio  'Pel  Peon- 
The  ohmic 


overpotentials 


Eq. 


(13) 

(5), 


Pohm>  Pohm-  Pohm-  Pohm  and  Pohm-  are  determined  by  the  elec- 
tronic  (ionic)  potential  differences  along  the  electronic  (ionic) 
current  flux  paths. 


2.1.2.  Effective  conductivities  and  effective  TPB  length 

The  effective  conductivity  of  the  /<-phase  material  (electronic  or 
ionic  conducting  material)  may  be  calculated  by  [37]: 


rreff  -  rrQ 

—  0  b 


Jk  -  fk 


1  +  $g/  (l  *^g)  ’/'k. 


(14) 


where  <pg  is  the  porosity  ratio,  ^  and  are  respectively  the  volume 
fraction  and  the  percolated  volume  fraction  threshold  of  the 
/<-phase  material.  Eq.  (14)  applies  only  for  \ pk  >  i/4-  °tff  *s  set  t0  zero 
when  fa  <  i/4-  i/4  is  calculated  as  [38], 


tfk/rk 


tfk/rk  +  (i  -  ik)  / rl 


=  1.764, 


(15) 


where  r2-  (z  =  k,  l )  is  the  radius  of  z-phase  particle,  Z  is  the  average 
coordination  number  of  all  solid  particles  and  set  to  6  [38,39]. 

The  intrinsic  conductivities  (<t£)  of  Ni,  LSM  and  YSZ  are  deter¬ 
mined  as  [7,40], 

<T°i  =  3.27  x  106  -  1065.37,  (16) 


4.2  x  107 


'LSM 


-exp 


-1150 


(17) 


<j?sz  =  6.25  x  104exp(  1°3°°j.  (18) 

The  effective  or  percolated  volume  specific  TPB  length  is 
determined  as: 

^TPB.eff  =  27rmin(red,re,)(sinP)n^dZed,e,PedPe,,  (19) 


where  n%  =  3(1  -  0g)^k/ 4rcr^  [7]  is  the  volumetric  number  density 
of  k  particles  in  a  composite  electrode,  6  is  the  angle  of  particle  contact 
(assumed  to  be  15°  [38]).  Zjy  is  the  average  number  of  /c-particles  in 
contact  with  a  /-particle  that  may  be  estimated  as  [33]: 


zk,l  = 


0.5(l  +rl/rf) 


ii/U 


E"  1  ik/rk 


(20) 


The  percolation  probability  of  the  electrode  (electrolyte)  parti¬ 
cles,  Ped  (Pei),  can  be  calculated  as  [38]: 


Pk 


^3.764  -Z^2yA 


(21) 


2.2.  Thermal  fluid  model 
2.2.1.  Mass  transport  equations 

This  model  describes  an  mtSOFC  unit  running  on  hydrogen  and 
air.  The  humidified  hydrogen  gas  is  supplied  as  fuel  in  the  anode. 
Air  is  supplied  as  oxidant,  consisting  of  two  components:  oxygen 
and  nitrogen. 

The  mass  transport  in  the  gas  channel  may  be  described  by  the 
mass  balance  equation: 

V ■  Nj  =  V •  (  -  DjVC,  +  q  u  )  =  Rj  (22) 

where  Nt  is  the  molar  flux  of  species  z,  Q  the  molar  concentration  of 
species  z,  u  the  vector  of  convection  velocity.  P2  is  the  reaction  rate 
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of  species  i  and  is  zero  in  the  gas  channel.  Du  the  diffusion  coeffi¬ 
cient  of  species  i,  is  the  same  as  the  binary  diffusion  coefficient  Dy  as 
the  fuel  and  air  each  consists  of  only  two  species.  The  binary  dif¬ 
fusion  coefficient  Dy  may  be  evaluated  as  [31  ], 


Dij  = 


3.198  x  lO-8!1  75 

,V3  ,  l/3\ 


P(v  i 


Vi  J 


1 

Mi 


1 


0.5 


(23) 


where  p  is  the  gas  pressure,  v  is  the  diffusion  volume  (6.12  x  10  6, 
13.1  x  1(T6, 16.3  x  1(T6  and  18.5  x  1(T6  m3  mol-1  for  H2,  H20,  02 
and  N2,  respectively  [41]),  M  is  the  molecular  mass  and  the  sub¬ 
script  i  {j )  refers  to  the  gas  species  i  (J). 

The  mass  transportjn  the  porous  electrodes  also  obeys  the  mass 
balance  equation,  V  •  Nf  =  Rim  The  rates  of  hydrogen,  oxygen  and 
steam  productions  are  calculated  respectively  as: 


B0,  Dff  and  Dfj?  are  calculated  by  the  following  equations 


[45,46]: 


Bn  = 


45(1  -*g) 
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(34) 


(35) 


where  e  is  the  porosity,  z  is  the  tortuosity  factor  and  rg  is  the  pore 
radius  determined  as  [33] 


le 

~2 F 

(25) 

ie 

_4F’ 
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ie 

z  2 r 
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where  ie  =  jTPB^TPB.eff  is  the  current  density.  The  molar  fluxes  of  gas 
species  in  the  porous  electrodes  should  be  described  by  the  dusty 
gas  model  that  takes  the  Knudsen  diffusion  into  account  [42].  For 
binary  flow,  the  molar  flux  may  be  written  as  [43]: 


=  Nf iffusion  - 

where 
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Here  x\  is  the  molar  fraction  of  species  i,Bo  is  the  permittivity. 
D?ff  and  Dfjf  are  respectively  the  effective  binary  diffusion  coeffi¬ 
cient  and  effective  Knudsen  diffusion  coefficient  of  species  i.  Ct0 ,t  is 
the  total  molar  concentration  of  species  in  the  electrode,  fi  is  the 
effective  viscosity  coefficient  of  the  mixture  fluid  that  may  be  cal¬ 
culated  as  [35]: 


ii  >7 1-vV 


(30) 


where  m  is  the  number  of  component  in  the  fluid  and  is 
a  dimensionless  number  determined  by: 


,  1  L  Mi 

iJ  ~  V8\  +Mj 


-0.5 


1+1^ 


0.5  , 

Mi  \  (  Mi 


0.25- 


(31) 


Here  pi  {pj)  is  the  viscosity  coefficient  of  species  i  (J),  which  may 
be  obtained  from  the  Sutherland  law  [44]  : 


Pi  =  Mref 


Tref  +  cref  (  T 


1.5 


(32) 


'  cref  \Jref  j 

where  pve f  is  the  reference  viscosity  at  the  reference  temperature, 
Tref,  and  cre f  is  the  Sutherland’s  constant. 


2  1 


1 


Y  _  _ _ _ 

3  1  —  ^Aed  Aed  +  ^Ael/^el 


(36) 


2.2.2.  Momentum  transport 

The  momentum  transport  in  the  fuel  channel  is  described  by  the 
Navier-Stokes  equation: 


T 

7-(fl(vu  +  (Vu)  )) 


-  Vp  =  P I  U  •  V  u , 


V-  u  =  0. 

Here  p  is  the  gas  density  given  by: 


(37) 

(38) 

(39) 


The  momentum  transport  in  the  porous  electrodes  is  modeled 
by  the  Brinkman  equation  that  accounts  for  the  viscous  transport  in 
the  momentum  balance  and  treats  both  the  pressure  and  the  flow 
velocity  vector  as  independent  variables  [47]: 


tu=-vP+v.(i(vu  +  (v„)r))-v-(; 


2/U, 

3^ 


V  ul  .  (40) 


2.2.3.  Heat  transport 

The  energy  conservation  equation  in  this  model  is  written  as: 


V-  - 


(-keffvr)  +  pCp  u 


vr  =  q, 


(41) 


where  Q  is  the  heat  source  that  includes  the  Ohmic  heat  source, 
Qohm,  activation  heat  source,  Qact,  and  the  entropy  heat  source,  Qentr- 
The  heat  sources  may  be  expressed  as: 


a  [2 

n  —  Al  _l  10 

xohm  —  _  “t  5 
(7ei  (T[0 

Qact  —  lelMact? 

Gemr  =  »el(-^)- 


(42) 

(43) 

(44) 


/<eff  in  Eq.  (41 )  is  the  effective  thermal  conductivity  of  the  porous 
material  and  determined  as: 


/<eff  —  0g/<f  +  ^1  0g^/<s, 


(45) 
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where  ks  is  the  thermal  conductivity  of  the  solid  material  and  kf  is 
the  effective  thermal  conductivity  of  the  fluid  that  is  calculated 
as  [41]: 


n 


kf  =  E 

i  =  1 


Xjk, 

ZjUxAj 


(46) 


where  Ay  is  evaluated  by  [41  ], 


(47) 


ki  in  Eq.  (46)  is  the  thermal  conductivity  of  species  i  and  deter¬ 
mined  as  [41]: 


Qrad  =  *rad(7?w-r4)>  (54) 

where  Tfw  is  the  temperature  of  the  furnace  wall.  erad  is  the  surface 
emissivity  of  the  cathode. 

Combining  both  convective  heat  transfer  and  radiation,  the  heat 
transport  at  the  interface  between  the  air  flow  and  the  mtSOFC 
cathode  is  described  by: 

_  n  ■  (  -  kWT  +  pCp  u  T)  =  h(Text  -  T)  +  £rad  (t^  -  T4) .  (55) 

To  be  precise,  the  material  properties  and  operating  parameters 
used  in  this  study  are  listed  in  Table  2. 

2.3.  Structural  mechanics  model 


^  =  bf  x  10 -2+bf  x  10  5T  +  bf  x  1CT8T2.  (48) 

The  required  parameters,  b  f ,  bf  and  bf,  are  obtained  by  fitting 
the  experimental  results  [48]  and  are  listed  in  Table  1. 

Cp  in  Eq.  (41)  is  the  effective  specific  heat  of  the  fluid  and 
evaluated  as  [41]: 

Cp  =  J>4,  (49) 


where  Clp  is  the  specific  heat  of  species  i  and  calculated  as  [41]: 

Cp  =  aj  +af  X  1(T3T  +  af  x  10 -6T2.  (50) 

The  required  parameters,  af ,  af  and  af,  are  obtained  by  fitting 
the  experimental  data  [48]  and  listed  in  Table  1. 

At  the  interface  between  the  air  flow  and  the  mtSOFC  cathode, 
the  heat  transport  due  to  the  convective  heat  transfer  and  radiation 
should  be  considered.  The  convective  heat  transfer  between  the  air 
and  the  cathode  outer  surface  may  be  described  with  the  Newton’s 
law  of  cooling: 

Q.con  =  h(Text-T),  (51) 


where  Tex t  is  the  external  temperature  of  air,  or  the  temperature  of 
the  furnace  wall  in  this  study,  h  is  the  heat  transfer  coefficient 
estimated  as  [29]: 


h  = 


Nukj 
d  ’ 


(52) 


where  d  is  the  hydraulic  diameter  and  Nu  is  the  Nusselt  number 
calculated  as  [48]: 


Nu  =  3.66 


0.0668  (d/l)RePr 
1+0.04  ((d/l)RePr) 


2/3  • 


(53) 


Here  Pr  is  the  Prandtl  number,  Re  is  the  Reynolds  number  and  l  is 
the  length  of  the  mtSOFC  tube. 

The  radiative  heat  transfer  between  the  cathode  outer  surface 
and  the  furnace  wall  may  be  calculated  as: 


This  study  focuses  on  both  the  residual  stress  determined  by  the 
material  manufacturing  process  and  material  properties  and  the 
operating  thermal  stress  based  on  the  temperature  profile  obtained 


Table  2 

Fuel  cell  material  properties  and  operating  parameters. 


Fuel  flow  rate  ( u0 ) 

1.06  ms-2  [26] 

Fuel  or  air  pressure  (Po) 

101,325  Pa  [26] 

Operating  temperature  (Top) 

1073  K  [26] 

Operating  voltage  (Vcen) 

0.7  V 

Molar  air  composition 

21%  02  and  79%  N2  [26] 

Molar  fuel  composition 

97%  H2  and  3%  H20  [26] 

Tortuosity  (t) 

3  [7] 

Thermal  conductivity  of  anode 

2  W  ITT2  s- 1  [49,50] 

Thermal  conductivity  of  YSZ 

2  W  irr2  s- 1  [49,50] 

Thermal  conductivity  of  LSM 

4  W  irr2  s- 1  [49,50] 

Activation  energies  for  the  anode  (EHl ) 

120  kj  mol-1  [36] 

Activation  energies  for  the  cathode  ( Eq2  ) 

130  kj  mol  1  [36] 

Forward  reaction  symmetric  factor  for  anode  (af) 

1.5  [7] 

Reverse  reaction  symmetric  factor  for  anode  (ft) 

0.5  [7] 

Forward  reaction  symmetric  factor  for  cathode  (af) 

0.65  [7] 

Forward  reaction  symmetric  factor  for  cathode  (ft) 

0.35  [7] 

Exchange  transfer  current  density  of  anode  (jg  ref ) 

8  x  10"3  Airr1  [7] 

Exchange  transfer  current  density  of  cathode  (j £  ref ) 

1.25  x  10-4  A  m_1  [7] 

Ni  radius  (rNi) 

0.33  x  lO"6  m 

YSZ  radius  (rY sz) 

0.3  x  lO" 6  m 

LSM  radius  (rLsM) 

0.3  x  lO"6  m 

Volume  fraction  of  Ni 

0.445  [26] 

Volume  fraction  of  LSM  (^lsm) 

0.479  [26] 

Porosity  of  anode  (<pg) 

0.4 

Porosity  of  cathode  (<pg) 

0.4 

Prandtl  number  for  air  (Pr) 

0.7  [51] 

Thermal  surface  emissivity  (erad) 

0.3  [52] 

Anode  length 

30  x  10"3  m  [26] 

Cathode  length 

11.8  x  10  3  m  [26] 

Electrolyte  length 

20  x  10"3  m  [26] 

Width  of  anode  current  collector 

5  x  10"3  m  [26] 

Anode  thickness 

0.15  x  10  3  m  [26] 

Cathode  thickness 

20  x  10-6m  [26] 

Electrolyte  thickness 

12  x  10-6m  [26] 

Tubular  interior  diameter 

0.85  x  lO'3  m  [26] 

Furnace  interior  diameter 

6  x  10"3  m  [26] 

Anode  and  electrolyte  stress  free  temperature 

1473  I<  [15] 

Cathode  stress  free  temperature 

1373  K  [70] 

Table  1 

Parameters  for  calculating  the  viscosities,  specific  heats  and  thermal  conductivities  of  gas  species  [48]. 


Gas 

For  viscosity 

For  specific  heat 

For  thermal  conductivity 

Prep'S) 

+ef 

Gef 

a } 

i 

a? 

i 

a3 

i 

b] 

bf 

bf 

h2 

8.41  x  10"6 

273  K 

96.7  K 

29.09 

0.836 

-0.327 

5.61 

39.36 

61.7 

h2o 

1.70  x  10"6 

273  K 

861  K 

30.0 

10.7 

-2.02 

-0.27 

6.64 

1.38 

02 

19.2  x  lO"6 

273  K 

139  K 

36.16 

0.845 

-0.749 

0.21 

8.31 

-1.41 

N2 

16.6  x  10"6 

273  K 

107  K 

27.32 

6.23 

-0.950 

1.77 

4.22 

0.506 
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from  the  above  described  thermal-fluid  electrochemical  model. 
Thermal  stress  was  calculated  using  the  solid  mechanics  model.  It  is 
assumed  that  all  the  materials  of  the  SOFC  undergo  elastic  defor¬ 
mation  when  subject  to  thermal  loads.  In  this  model,  no  shrinkage 
is  considered  during  the  reduction  of  NiO  into  Ni  and  the  co¬ 
efficients  of  thermal  expansion  (CTEs)  of  NiO/YSZ  and  Ni/YSZ  are 
considered  to  be  equal.  This  model  is  often  used  in  the  open  liter¬ 
ature  [53]  and  adopted  to  calculate  the  stress  distribution  of  the 
Ni/YSZ  anode  in  the  work. 

2.3.1.  The  stress-strain  relation 

In  the  solid  mechanics  model,  the  stress— strain  constitutive 
relation  may  be  written  as, 


a  =  D?el  +  a0, 


(56) 


Ccom  + 1 / (G1  - C2 )  +6\p2 (3/<2  +2G2)/ (5G2 (3/^ +4G2 ))  ’  (62) 

Here  JCCOm  and  GCOm  are  the  bulk  modulus  and  shear  modulus  of 
the  dense  composite  material,  respectively.  K 2  (/<i)  and  G2  (Gi)  are 
respectively  the  bulk  and  shear  modulus  of  the  matrix  (impurity)  in 
the  composite  material. 

The  effective  elastic  and  shear  modulus  for  the  porous  material 
are  given  by  the  approach  developed  by  Ramakrishnan  and  Aru- 
nachalam  [55],  effective  Young’s  modulus  Eeff  and  shear  modulus 
Geff  are  determined  as: 


(1-+)2 

£eff  “  £°1  +(2-3i>0)^g 


(63) 


where  a  is  the  stress  tensor,  D  is  the  elasticity  matrix,  fei  is  the 
elastic  strain,  oo  is  the  initial  stress  or  the  residual  stress.  In  this 
work,  the  gravitational  forces  are  assumed  to  have  a  negligible 
contribution  to  the  overall  stress  distribution  and  the  materials 
undergo  free  deformations. 

For  an  isotropic  material  in  axial  symmetry,  D  is  defined  as  [14]: 


D  = 


(1  +  u)(l  -2v) 


1  -  v  v  v  0 
v  1  -  v  v  0 
v  v  1  -  v  0 

0  0  0  1 


1  -2v 


(57) 


Geff  =  Q) 


1  +  (1 1  -  19u0)/(4  -  4v0)cp„ 


(64) 


where  the  subscript  “0”  refers  to  the  property  of  the  dense  material. 
Inserting  Eqs.  (63)  and  (64)  into  Eq.  (58),  the  effective  Poisson  ratio, 
ueff,  for  the  porous  material  may  be  calculated  as: 


ueff  =  7 


1  4vq  +  3(/)g  —  7\)q  (j)g 


4  1+  2  (j)^  —  3  Vq(J)q 


(65) 


The  effective  CTE  of  the  composite  electrode  may  be  calculated 
as  [56]: 


where  E  is  the  Young’s  modulus  and  v  is  the  Poisson’s  ratio  of  the 
material  defined  by: 

“-sr1’  (58) 

where  G  is  the  shear  modulus  of  the  material. 

The  elastic  strain,  fei>  is  a  component  of  the  total  strain,  ?,  that 
also  includes  the  thermal  strain,  ?th.  and  the  initial  strain,  f0, 


=  <*1^1  Kl  ±02^2^2 

faK-l  +fal<2 


(66) 


The  effect  of  porosity  on  CTE  is  neglected  as  studies  have  shown 
that  there  is  no  significant  influence  of  porosity  on  the  thermal 
expansion  [57,58].  The  mechanical  properties  of  the  dense  mate¬ 
rials  are  shown  in  Table  3. 


2.4.  Failure  probability  analysis 


?  —  ?el+?th  +  ?0-  (59) 

The  thermal  strain  may  be  calculated  as, 

?th  =  a(T-T{),  (60) 

where  a  is  for  the  CTE  of  the  material,  T  is  the  local  temperature 
obtained  from  the  thermal-fluid  model,  and  Tf  is  the  stress  free 
temperature  at  which  there  is  no  stress  accumulated  in  the 
material. 

Combining  Eq.  (59)  and  Eq.  (60)  and  considering  the  relation¬ 
ship  between  the  initial  strain  and  initial  stress,  (70  =  Df0,  Eq.  (56) 
may  be  rewritten  as, 

O  =  D(S  -  ?th  -  fo)  +  ffo  =  £>[£  -  a(r  -  7f)j .  (56a) 


The  mechanical  failure  in  a  ceramic  component  occurs  due  to 
the  brittle  nature  of  the  material.  Meanwhile,  the  ceramics  exhibit 
a  statistical  strength  distribution  [52].  Therefore,  analyzing  the 
failure  probability  of  the  ceramic  component  is  necessary.  The 
Weibull  approach  of  failure  analysis  is  adopted  here  to  estimate  the 
risk  flaws  of  the  ceramic  component.  The  Weibull  method  [65] 
calculates  the  failure  probability,  Pf,  of  the  ceramic  component 
with  the  magnitude  of  the  tensile  stress,  ow 

v 

where  00  is  the  Weibull  strength  at  the  reference  volume,  Vo,  V  is 
the  volume  of  the  ceramic  component  and  m  represents  the  Wei¬ 
bull  modulus  which  corresponds  to  a  shape  parameter.  For 


2.3.2.  Material  properties 

For  the  mtSOFC  studied  here,  YSZ  electrolyte  is  coated  on  porous 
Ni-YSZ  anode  and  LSM— YSZ  is  used  as  cathode.  Mechanical  prop¬ 
erties  for  the  dense  anode  and  cathode  are  calculated  using  the 
relations  given  for  heterogeneous  materials  based  on  the  “com¬ 
posite  sphere  method”  (CSM)  [54]: 

l/q 

Kcom  =  I<2  +  i/(Ki  _k2)  +  3^/(3/^  +4G2)’  (61 3 


Table  3 

Mechanical  properties  of  dense  materials. 


Material 

Young’  modulus 
(GPa) 

Poisson 

ratio 

CTE 

(10-eK-i) 

Density 
(103  kg  m"3) 

YSZ 

215  [13] 

0.317  [13] 

10.3  [59] 

6.04  [60] 

YSZ  (800  °C) 
Ni 

185  [13] 

200  [61] 

0.313  [14] 

16.9  [59] 

8.9  [14] 

Ni  (800  °C) 
LSM 

171  [62] 

95  [63] 

0.32  [14] 

12.4  [64] 

6.57  [60] 
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a  material  subjecting  to  a  multi-axial  stress,  the  total  survival 
probability  may  be  calculated  as  the  product  of  the  survival  prob¬ 
abilities  determined  by  the  three  principal  stresses,  crz.  As  a  result, 
the  total  failure  probability  is  calculated  as: 


pm 


-'-nL.«p(-/(S) 


mdV\ 

Vo/ 


(68) 


The  three  principal  stresses,  04,  02  and  (73,  are  the  components  of 
the  stress  tensor  when  the  basis  is  changed  in  such  a  way  that  the 
shear  stress  components  become  zero.  Moreover,  the  principal 
stresses  are  chosen  such  that  04  >  02  >  03.  The  stress  is  tensile  if 
G[  >  0,  and  the  stress  is  compressive  if  a\  <  0.  The  Weibull  analysis 
only  applies  for  a  tensile  stress.  0[  =  0  is  used  in  calculating  Eq.  (68) 
where  a  compressive  stress  (d/  <  0)  is  encountered  in  the  stress 
field.  The  Weibull  parameters  for  the  SOFC  materials  studied  here 
are  shown  in  Table  4  [13,53,66]. 

When  subjecting  to  a  compressive  stress,  the  material  survives  if 
the  stress  is  below  its  compressive  strength.  Otherwise,  the  mate¬ 
rial  fails  mechanically.  The  compressive  strengths  of  the  YSZ  elec¬ 
trolyte  and  the  LSM-YSZ  composite  cathode  are  listed  in  Table  5. 


3.  Numerical  method 


3.1.  Geometric  model 

Fig.  la  shows  an  anode-supported  mtSOFC.  The  inner  layer  is  the 
porous  anode  where  fuel  flows  through  the  cell  tube.  The  outer 
layer  is  the  porous  cathode  where  air  flows  around.  The  middle 
layer  is  the  dense  YSZ  electrolyte  that  conducts  oxygen  ions  but  not 
the  electrons.  There  is  an  axial  symmetry  in  the  mtSOFC  and  the 
numerical  model  is  two  dimensional  (2D)  in  nature.  The  2D  com¬ 
putational  domain  is  illustrated  in  Fig.  lb.  The  actual  three 
dimensional  geometry  of  the  mtSOFC  may  be  obtained  by  revolving 
the  computational  domain  around  the  symmetry  axis. 

3.2.  Boundary  conditions 

The  boundary  conditions  are  critically  important  for  the  solu¬ 
tions  of  the  above  described  partial  differential  equations  (PDEs)  for 
the  thermal  fluid  electrochemistry  model  and  the  structural  me¬ 
chanics  model  and  should  be  specified  for  all  relevant  boundaries. 
Table  6  summarizes  the  boundary  settings  for  the  multi-physics 
electrochemistry  model.  The  boundary  conditions  for  the  struc¬ 
tural  mechanics  are  set  to  be  free  deformation. 


3.3.  Numerical  solution 

The  PDEs  for  the  thermal  fluid  electrochemistry  model  and  the 
structural  mechanics  model  with  the  appropriate  boundary  con¬ 
ditions  were  solved  using  the  finite  element  commercial  software 
COMSOL  MULIPHSICS  Version  3.5a  [47].  The  COMSOL  stationary 


Table  4 

Weibull  parameters  for  SOFC  material  used  in  this  model  [13,54,67]. 


Material 

Weibull 
modulus,  m 

Weibull 
strength,  a0 
(MPa) 

Reference 
volume,  V0 
(10-9  m3) 

Anode  (298  I<) 

11.8 

187 

0.578 

Anode  (1073  K) 

6.1 

115 

1.0a 

Electrolyte  (298  K) 

7.0 

446 

0.35 

Electrolyte  (1073  I<) 

8.0 

282 

0.270 

LSM  (298  K) 

7.0 

52 

1.210 

LSM  (1073  K) 

4.0 

75 

2.810 

a  Assumed. 


Table  5 

Compressive  strength. 


Material 

Room 

Operating 

temperature 

temperature 

Electrolyte 

1  GPa  [67] 

235  MPa  [68] 

Cathode 

100  MPa  [69] 

23.5  MPa3 

a  Assumed. 


nonlinear  solver  uses  an  affine  invariant  form  of  the  damped 
Newton  method  to  solve  the  discretized  PDEs  with  a  relative  con¬ 
vergence  tolerance  of  1.0  x  10  6.  The  two  sets  of  PDEs  for  the  multi¬ 
physics  electrochemistry  model  and  the  structural  mechanics 
model  are  solved  sequentially:  the  thermal  fluid  electrochemistry 
model  is  solved  first  and  then  the  predicted  temperature  field  is 
used  in  the  solid  mechanics  model.  Unless  explicitly  stated  other¬ 
wise,  the  standard  parameter  sets  shown  in  Tables  1-6  are  used  as 
the  defaults  in  the  numerical  simulations  throughout  this  work. 

4.  Results  and  discussion 

4.1.  Performance  of  mtSOFC  with  the  default  material  parameters 

4.1.1.  Validation  of  the  numerical  model 

An  mtSOFC  with  Ni-YSZ  composite  anode,  YSZ  electrolyte  and 
LSM-YSZ  cathode  was  manufactured  and  its  I-V  relations  were 
measured  at  the  temperatures  of  600,  700  and  800  °C  [26].  In  the 
experiment,  the  current  was  collected  only  from  one  side  of  the 
anode  layer,  instead  of  from  both  sides  of  the  anode  tube  as  indi¬ 
cated  in  Fig.  1.  For  comparison,  the  numerical  model  also  used  the 
one-side  anode  current  collection  mode.  The  parameters  shown  in 
Tables  1-5  were  used.  Most  of  these  parameters  were  obtained  by 
previous  experimental  measurements  on  pSOFCs  [32]  or  deduced 


Fuel  outlet 

Fig.  1.  Schematic  of  mtSOFC  (a)  three-dimensional  structure,  (b)  two-dimensional 
computational  domain. 
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Table  6 

Boundary  conditions  for  the  thermal  fluid  electrochemistry  model. 
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Boundary 

Charge  transport 

Mass  transport 

Momentum  transport 

Heat  transport 

Fuel  channel  inlet 

H2  and  H20  concentration 

Laminar  flow  ( u0 ) 

Temperature  (Top) 

Fuel  channel  outlet 

Convective  flux 

Pressure,  no  viscous  stress  (P0) 

Convective  flux 

Interconnect— anode  interface 

Reference  potential  (. E0 ) 

Insulation 

No  slip 

Heat  flux  (Eq.  (55)) 

Air  channel— Cathode  interface 

Reference  potential  (Vceii) 

02  and  N2  concentration 

Heat  flux  (Eq.  (55)) 

Air  channel— electrolyte  interface 

Electric  insulation 

Insulation 

Heat  flux  (Eq.  (55)) 
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from  fitting  the  experiments  [7,32].  The  main  changes  concern  the 
particle  sizes  of  the  electrode  materials.  There  are  both  Ni-YSZ 
support  layer  and  functional  layer  in  the  measured  pSOFC  [32], 
while  there  is  only  one  single  Ni-YSZ  anode  layer  in  the  measured 
mtSOFC  [26].  Similarly,  there  are  two  cathode  layers  in  the  pSOFC 
and  one  cathode  layer  in  the  mtSOFC  [26,32].  Since  the  data  for  the 
particle  sizes  were  not  reported  in  the  mtSOFC  experiment  [26],  our 
numerical  simulation  of  the  mtSOFC  model  uses  the  particle  sizes 
(Table  2)  that  are  intermediate  to  the  values  used  for  the  pSOFC 
model  [7,31]  as  a  compromise.  Similarly,  intermediate  values  are 
used  for  the  electrode  porosities.  In  addition,  the  experiment  on 
pSOFC  found  large  contact  resistance  between  the  material  layers 
[32].  As  the  inter-layer  contact  resistance  was  not  examined  for 
mtSOFC,  the  contact  resistance  was  adjusted  to  fit  the  experimental 
mtSOFC  I—V  curves.  The  fit  gives  the  total  contact  resistance  of 
0.04  Q  cm2,  0.05  Q  cm2  and  0.07  Q  cm2  for  800  °C,  700  °C  and 
600  °C,  respectively.  These  values  are  smaller  than  the  corre¬ 
sponding  pSOFC  ones  [32]  and  are  quite  reasonable  as  there  are 
more  material  layers  in  the  pSOFC  than  in  the  mtSOFC. 

The  theoretical  and  experimental  I—V  relations  of  mtSOFC  at 
different  working  temperatures  are  compared  in  Fig.  2a.  Good 
agreement  between  the  theory  and  the  experiment  is  found,  dem¬ 
onstrating  that  the  present  thermal  fluid  electrochemistry  model  is 
capable  of  providing  quality  predictions.  It  also  indicates  that  the 
material  property  parameters  deduced  mostly  from  independent 
experimental  and  theoretical  results  are  of  general  applicability. 

4.1.2.  Effect  of  anode  current  collection  on  the  electrochemical 
performance 

The  electrochemical  performance  of  mtSOFC  shown  in  Fig.  2a  is 
much  poorer  than  its  pSOFC  counterpart.  One  major  contributing 
factor  is  the  current  collection  mode.  The  contour  of  the  current 
density  distribution  shown  in  Fig.  2b  gives  the  hint  for  the  poor 
performance.  All  current  has  to  pass  through  the  cross  section 
determined  by  the  anode  layer  thickness  (and  the  perimeter  of  the 
mtSOFC  tube)  and  the  average  length  of  the  current  passage  equals 
roughly  to  the  half  of  the  tube  (anode)  length.  If  the  current  is 
collected  by  both  sides  of  the  anode  current  collector,  the  cross 
section  area  will  be  increased  by  a  factor  of  two  and  the  length  for 
current  passage  will  be  reduced  by  half  roughly,  in  comparison  with 
the  case  of  one-side  current  collection.  The  ohmic  polarization  is 
thus  substantially  reduced  and  an  increased  cell  performance  is 
expected.  Indeed,  significantly  improved  cell  performance  (Fig.  3a) 
is  obtained  by  numerical  simulation  for  the  same  mtSOFC  with  the 
current  collection  by  both  sides  of  the  anode.  Comparison  of  Figs.  3a 
and  2a  show  that  the  overall  current  output  is  increased  by  about 
50%,  demonstrating  the  critical  importance  of  the  current  collection 
mode.  As  the  current  collection  by  both  outer  surfaces  of  the  anode 
may  be  easily  made  experimentally,  we  focus  on  this  current  col¬ 
lection  mode  in  the  following  numerical  simulations.  The  cell  with 
the  default  parameters  and  this  current  collection  mode  is  refer¬ 
enced  as  the  standard  cell  below. 


4.1.3.  Temperature  profile 

Fig.  4  shows  the  temperature  profile  of  mtSOFC  when  operating 
at  the  fuel  inlet  temperature  of  1073  K  and  the  cell  voltage  of  0.7  V. 
The  temperature  varies  from  1073  K  to  1090.2  K  for  the  cell  with 
a  tube  length  of  3  cm.  As  shown  in  Fig.  4,  the  two  ends  of  the  cell 
tube  have  a  lower  temperature  than  the  middle  of  the  cell  tube  and 
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Fig.  2.  Electrochemical  performance  of  mtSOFC  with  the  current  collected  by  one 
side  of  the  anode  (a)  comparison  of  the  theoretical  and  experimental  I—V  results; 
(b)  streamline  of  the  electronic  current  density  distribution  in  the  electrodes  for 
T  =  800  °C  and  Vcen  =  0.7  V.  The  thickness  direction  (r-direction)  is  magnified  for 
easy  display. 
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Fig.  3.  mtSOFC  with  the  current  collected  by  both  sides  of  the  anode  (a)  I-V  relations; 
(b)  the  electronic  current  density  distribution  for  T  =  800  °C  and  Vcen  =  0.7  V. 

the  hotspot  is  located  closer  to  the  fuel  exit  end  of  the  cell.  The 
temperature  variation  of  about  17  I<  along  the  cell  tube  of  3  cm  is 
rather  mild  in  comparison  with  that  found  for  pSOFC  [34].  The 
relatively  uniform  temperature  distribution  of  mtSOFC  indicates 
that  mtSOFC  may  operate  at  a  higher  temperature  for  the  increased 
power  output  than  its  pSOFC  counterpart. 

4.1.4.  Stress  distributions  for  T0  =  298  I<  and  Top  =  1073  I< 

The  residual  stress  in  the  cell  at  room  temperature  depends  on 
the  cell  fabrication  process.  For  consistency,  we  consider  the 
mtSOFC  made  by  Yang  et  al.  [26]  and  the  relevant  fabrication  pro¬ 
cesses  may  be  summarized  as  the  following.  The  electrolyte  layer 
was  coated  on  to  the  anode  surface  using  a  vacuum  assisted  coating 
technique  and  then  co-fired  at  1723  I<  for  10  h  to  form  a  thin  and 
dense  layer.  The  cathode  materials  were  coated  on  to  the  electrolyte 
by  brush  printing  and  sintered  at  1423  K  for  3  h.  With  these  fab¬ 
rication  processes,  the  zero  stress  temperatures  suggested  for  the 
anode,  electrolyte  and  cathode  are  1473  K,  1473  K  and  1373  K, 
respectively  [15,70].  We  assume  that  the  residual  stresses  are 
induced  in  the  anode/electrolyte/cathode  structures  only  due  to  the 
CTE  mismatch  between  different  layers.  The  thermal  stress  is 
introduced  by  the  CTE  mismatch  and  the  temperature  profile. 
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Fig.  4.  Temperature  profile  of  mtSOFC  operated  at  800  °C  and  0.7  V. 

Fig.  5a  shows  the  residual  stress  distribution  of  the  cell  at  room 
temperature.  The  blue,  red  and  green  lines  stand  for  the  first,  sec¬ 
ond  and  third  principal  stresses,  respectively.  As  shown  in  Fig.  5a, 
the  stress  in  the  anode  is  mostly  tensile,  while  the  stresses  in  the 
cathode  and  the  electrolyte  are  mostly  compressive.  As  the  cell  is 
manufactured  at  high  temperature,  it  suffers  a  great  residual  stress 
at  room  temperature.  The  maximum  residual  compressive  stress  is 
375  MPa  for  the  electrolyte  and  26  MPa  for  the  cathode.  The  re¬ 
sidual  tensile  stress  for  the  anode  is  around  35  MPa. 

Fig.  5b  shows  the  stress  distribution  in  the  middle  cross  section 
of  the  cell  at  the  operating  temperature  of  1073  K  and  the  operating 
voltage  of  0.7  V.  The  stress  distribution  is  calculated  with  the 
temperature  profile  of  the  cell  at  the  operating  condition.  Similar  to 
the  case  for  the  room  temperature,  the  stresses  in  the  operating  cell 
are  mostly  tensile  for  the  anode  and  mostly  compressive  for  the 
cathode  and  the  electrolyte.  The  magnitudes  of  the  thermal  stresses 
in  the  operating  cell  are  found  to  be  smaller  than  that  of  the  room 
temperature  residual  stresses.  This  is  due  to  that  the  stress  free 
temperatures  are  closer  to  the  operating  temperature  than  to  the 
room  temperature.  The  maximum  compressive  stress  is  about 
123  MPa  for  the  electrolyte  and  7  MPa  for  the  cathode.  The  max¬ 
imum  tensile  stress  of  the  anode  is  14  MPa.  Even  though  the 
thermal  stress  is  smaller  at  the  operating  temperature  than  at  room 
temperature,  the  mechanical  strength  of  the  electrode  is  also  lower 
at  the  operating  temperature.  As  to  be  seen  below,  the  failure 
probability  turns  out  to  be  higher  at  the  operating  temperature. 

4.2.  Effect  of  electrode  composition  on  electrochemical  performance 

As  seen  in  Section  2.1,  the  material  compositions  such  as  the 
volume  fractions  and  the  particle  sizes  of  Ni/LSM/YSZ  affect  both 
the  TPB  lengths  and  electrical  conductivities  of  the  electrodes. 
Consequently,  the  material  compositions  affect  both  the  activation 
and  ohmic  polarizations.  Here  we  apply  the  mathematical  model  to 
examine  the  effects  of  the  electrode  compositions  and  their  designs 
for  improving  the  electrochemical  performance. 

4.2.1.  Effect  of  electrode  composition  on  the  TPB  length 

As  shown  in  Eq.  (19),  the  anode  (cathode)  TPB  length  depends 
on  the  percolation  probability  of  Ni  (LSM)  and  YSZ,  the  average 
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Fig.  5.  Distributions  of  the  first,  second  and  third  principal  stresses  of  mtSOFC  along 
the  radial  direction  (a)  T  =  298  K,  (b)  Top  =  1073  K. 


coordination  number  between  Ni  (LSM)  and  YSZ,  and  the  number 
of  Ni  (LSM)  particles  per  unit  volume.  The  dependence  of  the  per¬ 
colated  TPB  length  on  the  volume  fraction  and  particle  size  of  Ni 
(LSM)  in  the  anode  is  shown  in  Fig.  6.  The  volume  fraction  of  Ni 
(LSM)  is  varied  from  0.3  to  0.7.  The  radius  of  the  Ni  (LSM)  particle 
changes  from  0.18  x  10-6  m  (0.21  x  10  6  m)  to  0.42  x  10-6  m 
(0.42  x  10  6  m),  while  the  radius  of  the  YSZ  particle  is  kept  at 
0.3  x  10  6  m.  The  maximum  lengths  of  the  anode  TPB  for  the  Ni 
radii  of  0.18  x  10-6  m,  0.33  x  10~6  m  and  0.42  x  10-6  m  are 
respectively  5.12  x  1012  m~2,  3.12  x  1012  m-2  and  2.26  x  1012  m  2, 
corresponding  to  the  Ni  volume  fractions  of  0.39,  0.53  and  0.59, 
respectively.  The  anode  TPB  length  normally  increases  with  the 
decrease  of  the  Ni  radius  due  to  the  increased  number  of  Ni  par¬ 
ticles  per  unit  volume.  However,  the  change  of  the  Ni  radius  also 
affects  the  coordination  numbers  (Eq.  (20))  and  the  percolation 
probabilities  (Eq.  (21)).  The  percolation  regime  shifts  to  smaller 
volume  fraction  of  Ni  when  the  Ni  particle  size  decreases.  Con¬ 
sequently,  the  Ni  volume  fraction  corresponding  to  the  maximal 


Fig.  6.  Variation  of  the  anode  TPB  length  with  the  Ni  volume  fraction  and  the  Ni  radius. 
The  radii  of  the  YSZ  particles  are  kept  at  0.3  x  10-6  m. 


TPB  length  decreases  with  the  decreased  Ni  radius,  as  observed.  The 
trend  is  the  same  for  the  cathode.  The  maximum  lengths  of 
the  cathode  TPB  for  the  LSM  radii  of  0.21  x  10-6  m,  0.3  x  10~6  m 
and  0.42  x  10^6  m  are  respectively  4.47  x  1012  m-2, 3.59  x  1012  m~2 
and  2.26  x  1012  m~2,  corresponding  to  the  LSM  volume  fractions  of 
0.4,  0.5  and  0.59,  respectively. 

4.2.2.  Effects  of  electrode  composition  on  electrical  conductivity 

The  effects  of  the  Ni  (LSM)  volume  fraction  and  particle  size  on 

the  electronic  and  ionic  conductivities  of  the  anode  (cathode)  as 
determined  by  Eq.  (14)  are  shown  in  Fig.  7  for  the  operating  tem¬ 
perature  of  1073  K.  The  electronic  conductivity  of  the  anode 
(cathode)  increases  with  the  increased  Ni  (LSM)  volume  fraction  or 
the  reduced  Ni  (LSM)  particle  size  due  to  the  increased  relative 
abundance  of  the  electronic  conducting  Ni  (LSM)  particles.  Con¬ 
versely,  the  ionic  conductivity  of  the  anode  (cathode)  decreases 
with  the  increased  Ni  (LSM)  volume  fraction  or  the  reduced  Ni 
(LSM)  particle  size  for  the  same  reason.  An  increased  electronic 
conductivity  is  helpful  for  reducing  the  ohmic  polarization  of  car¬ 
rying  the  electronic  current,  while  an  increased  ionic  conductivity 
reduces  the  ohmic  polarization  of  carrying  the  ionic  current, 
allowing  for  the  electrochemical  reaction  to  take  place  in  an 
increased  distance  from  the  electrolyte  surface.  The  two  competing 
factors,  together  with  the  consideration  about  the  TPB  length, 
should  be  properly  balanced  to  obtain  the  desirable  cell  perfor¬ 
mance,  as  to  be  discussed  below. 

4.2.3.  Effect  of  anode  composition  on  electrochemical  performance 

Fig.  8a  shows  the  average  current  densities  for  the  cell  operating 

voltages  of  0.6, 0.7  and  0.85  V  as  functions  of  the  Ni  volume  fraction. 
The  Ni  volume  fractions  for  producing  the  maximum  current 
densities  for  the  three  operating  voltages  are  all  about  0.65.  That  is, 
the  optimal  Ni  volume  fraction  for  the  maximal  current  output  is 
insensitive  to  the  operating  voltage.  However,  the  current  density 
changes  more  rapidly  with  the  Ni  content  for  the  lower  operating 
voltage  than  for  the  higher  operating  voltage.  For  example,  the 
current  densities  for  the  cell  voltage  of  0.85  V,  0.7  V  and  0.6  V  are 
respectively  93%,  87%  and  84%  of  their  maximal  when  the  Ni  vol¬ 
ume  fraction  is  changed  to  0.5  from  its  optimal  value  of  0.65.  The 
stronger  influence  of  the  Ni  volume  fraction  on  the  current  output 
of  the  lower  operating  voltage  is  understandable  as  the  current 
produced  by  the  lower  operating  voltage  is  larger  and  is  therefore 
more  sensitive  to  the  effect  of  the  Ni  content  on  the  electrical 
conductivity  and  the  TPB  length. 

Even  though  the  anode  electronic  conductivity  is  much  larger 
than  the  anode  ionic  conductivity  for  most  of  the  Ni  content 
examined  (Fig.  7),  the  optimal  electrode  composition  corresponds 
to  a  content  of  Ni  higher  than  that  of  YSZ.  This  indicates  that  the 
anode  electronic  conductivity  is  more  important  than  the  anode 
ionic  conductivity  for  the  current  output.  The  reason  is  that  the 
electrochemical  reaction  happens  at  the  TPBs  near  the  electrode- 
electrolyte  interface  and  the  electronic  current  path  is  much  lon¬ 
ger  than  that  of  the  ionic  current  path. 

The  combination  effect  of  the  volume  fraction  and  particle  size 
of  Ni  is  shown  in  Fig.  8b.  The  maximal  current  output  for  the 
Ni  radius  of  0.18,  0.33  and  0.42  x  10  6  m  is  obtained  with  the  Ni 
volume  fraction  of  0.55,  0.6  and  0.7,  respectively.  As  shown  in 
Fig.  8b,  the  maximal  currents  are  not  very  different  for  different  Ni 
radii.  For  example,  the  maximum  current  densities  for  the  Ni  radii 
of  0.42  X  10"6  m.  0.33  x  10~6  m  and  0.18  x  10"6  m  are 
10,553  A  nr2,  10,846  A  m“2  and  10,992  A  m“2,  respectively. 
However,  the  current  output  decreases  more  slowly  for  smaller  Ni 
particles  than  for  larger  Ni  particles  when  the  Ni  volume  fractions 
deviate  away  from  their  optimal  values.  The  current  density  for  the 
Ni  radius  of  0.18  x  10-6  m  is  about  88%  of  its  maximum  when  the  Ni 
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Fig.  7.  Influence  of  the  material  composition  on  the  electrode  electrical  conductivity  at  T  =  1073  K.  (a)  Electronic  conductivity  of  the  anode  vs.  the  Ni  volume  fraction  and  particle 
size;  (b)  ionic  conductivity  of  the  anode  vs.  the  Ni  volume  fraction  and  particle  size;  (c)  electronic  conductivity  of  the  cathode  vs.  the  LSM  volume  fraction  and  particle  size;  (d)  ionic 
conductivity  of  the  cathode  vs.  the  LSM  volume  fraction  and  particle  size. 


volume  fraction  is  0.4,  or  a  deviation  of  about  12%  from  the  optimal 
Ni  volume  fraction  of  0.55.  For  a  similar  amount  of  deviation,  the 
current  outputs  for  the  Ni  radii  of  0.33  x  10-6  m  and  0.42  x  10-6  m 
are  only  about  62%  and  16%  of  their  maxima,  respectively. 


Fig.  8.  Effect  of  the  anode  composition  on  the  cell  output  (a)  current  density  vs.  the  Ni 
volume  fraction  and  the  operating  voltage;  (b)  current  density  vs.  the  Ni  volume 
fraction  and  particle  size. 


Therefore,  using  small  Ni  particle  sizes  allows  for  a  large  adjustable 
range  of  the  Ni  content  while  maintaining  the  desirable  current 
production. 

4.2.4.  Effect  of  cathode  composition  on  electrochemical 
performance 

The  dependence  of  the  cell  current  output  on  the  operating 
voltage  and  the  LSM  volume  fraction  is  shown  in  Fig.  9a.  The  vol¬ 
ume  fraction  of  LSM  for  producing  the  maximum  output  is  about 
0.4.  Unlike  the  case  for  the  anode,  the  cathode  should  have  more 
YSZ  than  LSM  to  be  optimal.  That  is,  the  cathode  ionic  conductivity 
is  often  more  important  than  the  cathode  electronic  conductivity 
for  the  current  output.  This  is  characteristically  different  from  that 
of  the  anode.  The  phenomena  is  attributed  to  that  the  electro¬ 
chemical  reaction  in  the  cathode  is  much  slower  than  that  in  the 
anode.  High  cathode  ionic  conductivity  is  important  for  increasing 
the  reaction  region  and  reducing  the  activation  polarization. 

Above  the  optimal  LSM  volume  fraction  of  0.4,  the  current 
density  increases  with  the  decrease  of  the  LSM  volume  fraction. 
Below  the  optimal  LSM  content,  the  current  density  initially  de¬ 
creases  slowly  with  the  decrease  of  the  LSM  volume  fraction,  but 
drops  rapidly  when  the  LSM  volume  fraction  is  approaching  0.3. 
This  is  because  that  the  percolation  probability  of  the  LSM  particle 
network  is  low  for  the  LSM  volume  fraction  below  0.35,  resulting 
in  an  unacceptably  low  electronic  conductivity  of  the  cathode 
(Fig.  7c,  rLSM  =  0.3  pm). 

Fig.  9b  shows  the  combined  effect  of  the  volume  fraction  and 
particle  size  of  LSM  on  the  current  output.  The  maximal  current 
outputs  for  the  LSM  radii  of  0.21, 0.3  and  0.42  x  10-6  m  are  obtained 
with  the  LSM  volume  fraction  of  0.3,  0.4  and  0.5,  respectively. 
Compared  to  the  case  for  the  anode  Ni  composition,  the  depend¬ 
ence  of  the  maximal  current  density  on  the  LSM  particle  size  is 
relatively  large  due  to  the  relatively  strong  effect  of  the  TPB  on  the 
cathode  electrochemical  reaction.  The  maximum  current  density 
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Fig.  9.  Effect  of  the  cathode  composition  on  the  cell  output  (a)  current  density  vs.  the 
LSM  volume  fraction  and  the  operating  voltage;  (b)  current  density  vs.  the  LSM  volume 
fraction  and  particle  size. 


increases  from  7565  A  m_2-8498  A  m”2,  or  an  increase  of  12%, 
when  the  LSM  radius  is  decreased  from  0.42  x  10  6  m  to 
0.21  x  10~6  m.  It  is  desirable  to  use  small  LSM  particle  to  improve 
the  cell  output. 


4.3.  Effect  of  electrode  composition  on  thermal  mechanical 
performance 

4.3 A.  Effect  of  electrode  composition  on  thermal  mechanical 
properties 

The  compositions  affect  the  thermal  mechanical  properties  of 
the  electrodes  according  to  the  theory  described  in  Section  2.3.  The 
Poisson  ratios  of  pure  Ni,  YSZ  and  LSM  as  well  as  the  Young’s  module 
of  LSM  at  the  mtSOFC  operating  temperature  are  assumed  here  to 
be  the  same  as  the  ones  at  room  temperature.  Fig.  10a  (Fig.  10b) 
shows  the  variations  of  the  anode  (cathode)  Young’s  module  and 
Poisson  ratio  with  the  volume  fraction  of  Ni  (LSM).  As  the  Poisson 
ratios  of  Ni  (LSM)  and  YSZ  are  very  similar,  the  Poisson  ratio  of  the 
Ni-YSZ  composite  anode  (LSM-YSZ  composite  cathode)  remains 
almost  constant  with  the  variation  of  the  Ni  (LSM)  volume  fraction. 
In  the  other  hand,  the  Young’s  module  shows  a  larger  variation  with 
the  LSM  volume  fraction  for  the  cathode  than  with  the  Ni  volume 
fraction  for  the  anode  as  the  difference  in  the  Young’s  module  is 
larger  between  YSZ  and  LSM  than  between  YSZ  and  Ni. 

Fig.  10c  shows  the  anode  CTE  as  a  function  of  the  Ni  volume 
fraction.  As  the  CTEs  of  Ni  and  YSZ  are  very  different,  applying  Eq. 
(66)  to  the  Ni-YSZ  composite  anode  is  found  to  produce  sizable 
deviations  from  the  available  experimental  data  [56].  Numerical 
tests  show  that  the  mechanical  failure  probabilities  of  the  elec¬ 
trodes  are  sensitive  to  the  anode  CTE.  Therefore,  it  is  necessary  to 
improve  the  numerical  quality  of  the  anode  CTE  for  the  interested 
Ni  content  (i/'Ni  e  [0.3, 0.7]).  For  that  purpose,  the  CTE  of  the  Ni-YSZ 
composite  anode  for  the  structural  mechanics  model  is  calculated 
by  the  following  empirical  fitting  expression: 

CTEanode  =  CTEysz  +  (CTENi  “  CTEYSz),//Ni(0 +°-9,//Ni)  (66a) 

As  seen  in  Fig.  10c,  our  fitting  expression  of  Eq.  (66a)  agrees 
reasonably  well  with  the  experimental  result. 


Fig.  10.  Dependences  of  the  thermal  mechanical  properties  on  the  electrode  compositions  (a)  Young’s  module  and  Poisson  ratio  of  the  anode  vs.  the  Ni  volume  fraction;  (b)  Young’s 
module  and  Poisson  ratio  of  the  cathode  vs.  the  LSM  volume  fraction;  (c)  anode  CTE  as  a  function  of  the  Ni  volume  fraction  (the  experimental  data  are  taken  from  Ref.  [59]);  (d) 
cathode  CTE  vs.  the  LSM  volume  fraction. 
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The  cathode  CTE  as  a  function  of  the  LSM  volume  fraction  as 
determined  by  Eq.  (66)  is  shown  in  Fig.  lOd.  As  the  CTEs  of  LSM  and 
YSZ  are  only  moderately  different,  Eq.  (66)  should  provide  a  good 
estimate  of  the  CTE  of  the  LSM-YSZ  composite  cathode.  In  fact,  as 
shown  in  Fig.  lOd,  the  results  by  Eq.  (66)  for  any  LSM  volume 
fractions  are  in  good  agreement  with  the  results  by  the  analog  of 
Eq.  (66a)  for  the  cathode,  CTEcathode  =  CTEysz  +  (CTElsm  —  CTEysz) 
’//lsm(0.1  +  0.9^lsm)*  This  comparison  may  be  seen  as  indirect  evi¬ 
dence  that  the  cathode  CTE  is  determined  with  good  accuracy. 
Moreover,  Fig.  10c  and  d  seems  to  imply  that  Eq.  (66)  is  only  valid 
for  the  composite  consisting  of  two  components  with  similar  CTEs, 
while  Eq.  (66a)  is  valid  even  when  the  CTEs  of  the  two  components 
are  substantially  different.  However,  as  numerical  results  show  that 
the  mechanical  failure  probabilities  of  the  electrodes  are  not  sen¬ 
sitive  to  the  cathode  CTE  and  to  conform  to  the  mainstream  liter¬ 
ature,  Eq.  (66)  is  used  to  determine  the  CTE  of  the  composite 
cathode  in  the  following  study. 

4.3.2.  Effect  of  anode  composition  on  mechanical  failure 

As  described  in  Section  2.4,  the  tensile  failure  of  the  ceramic 
components  of  the  cell  is  analyzed  using  the  Weibull  statistical 
method,  while  the  mechanical  failure  occurs  when  the  maximum 
compressive  stress  exceeds  the  compress  strength  of  the  material. 
The  sizes  of  the  Ni  and  LSM  particles  do  not  affect  the  mechanical 
properties  in  the  model  used  here.  As  described  above,  the  particle 
sizes  do  affect  the  electrochemical  performance.  However,  their 
effect  on  the  temperature  profile  is  rather  moderate  due  to  the 
small  dimension  of  the  mtSOFC.  Therefore,  the  analysis  of  the 
failure  probability  is  only  slightly  affected  by  the  particle  sizes.  The 
following  analysis  uses  a  Ni  radius  of  0.18  pm  and  a  LSM  radius  of 
0.3  pm  as  the  representing  case. 

Fig.  11a  shows  the  tensile  failure  probabilities  of  the  electrodes 
at  room  temperature  as  functions  of  the  Ni  content  in  the  anode. 
The  failure  probability  of  the  anode  is  only  3.8  x  10-11  for  the  Ni 
volume  fraction  of  0.3,  but  increases  to  0.0019  when  the  Ni  volume 


fraction  is  increased  to  0.7.  It  has  been  proposed  that  the  failure 
probability  should  be  under  10-5  to  ensure  mechanical  stability 
[71  ].  Based  on  this  criterion,  the  Ni  volume  fraction  should  be  below 
0.5.  The  failure  probability  of  the  cathode  is  smaller  than  1CT9  for  all 
interested  Ni  content  and  the  Ni  content  is  inconsequential  to  the 
mechanical  stability  of  the  cathode. 

Fig.  lib  shows  the  compressive  stress  strength  ratios  of  the 
electrolyte  and  the  cathode,  the  ratios  of  the  maximum 
compressive  stresses  of  the  electrolyte  and  the  cathode  to  their 
compressive  strengths,  at  room  temperature.  For  the  range  of  the  Ni 
volume  fraction  examined,  the  ratios  are  smaller  than  1  and  the 
mechanical  fracture  is  not  expected.  Therefore,  all  interested  Ni 
volume  fraction  meets  the  requirement  for  the  mechanical  stability 
of  the  electrolyte  and  cathode. 

Fig.  12a  shows  the  failure  probabilities  of  the  electrodes  at  the 
operating  temperature  as  functions  of  the  Ni  volume  fraction.  The 
anode  failure  probability  increases  from  2.93  x  10~7  to  0.003  when 
the  Ni  volume  fraction  increases  from  0.3  to  0.7,  while  the  cathode 
failure  probability  remains  below  10-9.  The  numerical  result  on  the 
failure  probability  clearly  shows  that  increasing  the  Ni  content  may 
easily  bring  about  mechanical  flaw  in  the  anode,  agreeing  well  with 
experimental  observation  [27].  Considering  the  criterion  of  the 
failure  probability  of  10-5  for  the  mechanical  safety,  the  Ni  volume 
fraction  should  not  exceed  0.42. 

The  compressive  stress  strength  ratios  of  the  electrolyte  and  the 
cathode  at  the  operating  temperature  are  shown  in  Fig.  12b.  Unlike 
the  case  for  the  room  temperature,  the  mechanical  failure  is 
expected  for  the  electrolyte  when  the  Ni  volume  fraction  exceeds 
0.65.  As  the  mechanical  stability  of  the  anode  requires  a  lower  Ni 
volume  fraction,  the  Ni  content  is  determined  by  the  anode 
stability. 

4.3.3.  Effect  of  cathode  composition  on  mechanical  failure 

The  effect  of  the  LSM  content  on  the  failure  probabilities  of  the 
electrodes  and  the  compressive  stress  strength  ratios  of  the 


Fig.  11.  Dependence  of  the  mechanical  failure  at  room  temperature  on  the  anode 
composition  (a)  failure  probabilities  of  the  electrodes  vs.  the  Ni  volume  fraction; 
(b)  ratios  of  the  maximum  compressive  stresses  of  the  electrolyte  and  the  cathode 
to  their  compressive  strengths  as  functions  of  the  Ni  volume  fraction. 


Fig.  12.  Dependence  of  the  mechanical  failure  at  the  operating  temperature  on  the 
anode  composition  (a)  failure  probabilities  of  the  electrodes  vs.  the  Ni  volume  fraction; 
(b)  compressive  stress  strength  ratios  of  the  electrolyte  and  the  cathode  vs.  the  Ni 
volume  fraction. 
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Fig.  13.  Dependence  of  the  mechanical  failure  at  room  temperature  on  the  cathode 
composition  (a)  failure  probabilities  of  the  electrodes  vs.  the  LSM  volume  fraction; 
(b)  compressive  stress  strength  ratios  of  the  electrolyte  and  the  cathode  vs.  the  LSM 
volume  fraction. 

electrolyte  and  the  cathode  at  room  temperature  is  shown  in  Fig.  13. 
The  failure  probabilities  of  the  electrodes  decrease  with  the 
increase  of  the  LSM  content  due  to  that  the  increased  CTE  of  the 
cathode  matches  better  with  that  of  the  anode,  reducing  the  tensile 
stresses  of  both  the  electrodes  as  well  as  the  compressive  stress  of 
the  cathode.  The  electrode  failure  probability  is  smaller  than 
6  x  10  8  for  all  interested  LSM  volume  fraction  and  is  well  under 
the  failure  probability  of  10-5  for  mechanical  stability.  The  com¬ 
pressive  stress  strength  ratios  of  the  electrolyte  and  the  cathode  are 
less  than  0.4  for  all  interested  LSM  volume  fraction  and  are  safe 
mechanically. 

The  effect  of  the  LSM  content  on  the  failure  probabilities  of  the 
electrodes  and  the  compressive  stress  strength  ratios  of  the  elec¬ 
trolyte  and  the  cathode  at  the  operating  temperature  is  shown  in 
Fig.  14.  The  failure  probability  of  the  anode  changes  slowly  with  the 
LSM  content.  It  decreases  from  1.26  x  10-5  to  9.0  x  10-6  when  the 
LSM  volume  fraction  increases  from  0.3  to  0.7.  These  failure  prob¬ 
abilities  are  at  the  borderline  of  10-5  for  mechanical  security  and 
should  be  acceptable  in  practice.  The  compressive  stress  strength 
ratios  of  the  electrolyte  and  the  cathode  are  less  than  0.6  for  all 
interested  LSM  volume  fraction  and  are  safe  mechanically. 

As  seen  from  the  above  discussion,  the  LSM  content  is  incon¬ 
sequential  from  the  mechanical  stability  point  of  view.  The  LSM 
content  should  be  determined  by  the  optimal  electrochemical 
performance  of  the  cell. 

4.4.  Primitive  comparison  of  the  current  outputs  ofmtSOFC  and 
pSOFC 

mtSOFCs  have  long  been  recognized  for  their  high  performance 
on  the  thermal  shock  resistance,  fast  startup  and  thermal  cycling. 
However,  the  current  density  output  of  mtSOFC  is  widely  shown  to 
be  comparable  with  or  moderately  superior  over  its  tSOFC  coun¬ 
terpart,  but  is  substantially  inferior  to  the  pSOFC  [2,26].  Never¬ 
theless,  as  discussed  above,  the  current  output  of  mtSOFC  may  be 
significantly  increased  by  changing  the  current  collection  from  one 


Fig.  14.  Dependence  of  the  mechanical  failure  at  the  operating  temperature  on  the 
cathode  composition  (a)  failure  probabilities  of  the  electrodes  vs.  the  LSM  volume 
fraction;  (b)  compressive  stress  strength  ratios  of  the  electrolyte  and  the  cathode  vs. 
the  LSM  volume  fraction. 

side  of  the  anode,  which  is  the  main  method  currently  in  use  [2,26], 
to  both  sides  of  the  anode.  In  fact,  the  current  production  shown  in 
Fig.  3a  is  not  very  far  away  from  the  state-of  the-art  performance  of 
pSOFC  [32].  It  would  be  interesting  to  make  a  comparison  of  the 
current  outputs  of  mtSOFC  and  pSOFC  on  a  reasonably  equal  basis. 
However,  due  to  the  involvement  of  a  large  number  of  different 
design  parameters,  a  quality  comparison  of  the  mtSOFC  and  pSOFC 
performances  requires  a  specifically  devoted  effort  and  only 
a  primitive  comparison  is  attempted  here. 

In  Ref.  [32],  the  electrochemical  performance  of  pSOFC  is  sys¬ 
tematically  improved  over  that  of  a  standard  cell  by  optimizing  the 
electrolyte  thickness,  cathode  interlayer  thickness,  anode  support 
layer  thickness  and  anode  support  porosity,  while  the  performance 
of  their  standard  cell  is  already  very  close  to  the  best  reported  by 
other  pSOFC  groups  using  the  same  set  of  materials.  Therefore,  it  is 
reasonable  to  use  the  standard  cell  performance  of  Ref.  [32]  as  the 
benchmark  for  comparison  here.  The  standard  cell  of  mtSOFC 
defined  above  is  modified  in  order  to  compare  with  the  standard 
cell  of  pSOFC  in  Ref.  [32]  on  an  equal  basis.  The  changes  made  to  the 
standard  cell  of  mtSOFC  are:  8  pm  for  the  electrolyte  thickness,  26% 
for  the  cathode  porosity,  23%  for  the  anode  porosity.  These  data  are 
the  corresponding  parameters  for  the  standard  pSOFC  cell  and  the 
modified  standard  mtSOFC  cell  may  then  mimic  the  standard 
pSOFC  cell  better.  The  I-V  curves  of  the  mtSOFC  and  pSOFC  cells  are 
compared  in  Fig.  15. 

As  shown  in  Fig.  15,  the  power  densities  of  the  standard  mtSOFC 
and  pSOFC  cells  are  quite  similar  for  the  operating  temperature  of 
700  °C.  The  maximal  power  density  of  the  pSOFC  cell  is  higher  than 
that  of  the  mtSOFC  cell  for  T  =  800  °C,  but  is  lower  than  that  of 
mtSOFC  for  T  =  600  °C.  For  practical  applications,  SOFC  has  to  work 
at  a  relatively  high  voltage  (typically  70-80%  of  the  open  voltage)  in 
order  to  have  the  desirable  energy  efficiency.  In  such  cases,  the 
performances  of  the  mtSOFC  and  pSOFC  cells  are  essentially  the 
same,  as  may  be  easily  seen  in  Fig.  15.  Therefore,  it  is  concluded  that 
the  output  current  density  of  mtSOFC  with  proper  current  collec¬ 
tion  is  similar  to  instead  of  much  smaller  than  that  of  pSOFC. 
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Current  Density  (Acrrf2) 


Fig.  15.  Comparison  of  the  electrochemical  performances  of  the  standard  mtSOFC  and 
pSOFC  cells. 


recommended  to  obtain  both  high  electrochemical  performance 
and  mechanical  stability.  Regarding  the  cathode  material,  the  best 
electrochemical  performance  of  the  cell  is  often  obtained  with 
a  LSM  volume  fraction  of  30—50%.  The  optimal  LSM  content  de¬ 
creases,  while  the  overall  electrochemical  performance  increases, 
with  the  reduced  LSM  particle  size.  Increasing  the  LSM  content 
reduces  the  compressive  stress  of  the  cathode,  but  its  effects  on  the 
anode  and  electrolyte  stresses  are  very  limited.  Overall,  the  LSM 
content  is  inconsequential  on  the  mechanical  stability  and 
should  be  optimized  based  on  achieving  high  electrochemical 
performance. 

A  primitive  comparison  of  the  I—V  curves  of  mtSOFC  and  pSOFC 
is  made.  It  is  concluded  that  the  current  outputs  of  pSOFC  and 
mtSOFC  are  in  fact  comparable  if  a  proper  current  collection 
method  is  used  in  mtSOFC.  Therefore,  mtSOFC  is  a  promising 
technology  with  both  the  benefits  of  pSOFC  for  high  current  density 
and  tSOFC  for  thermal  cycling  endurance. 


5.  Summary 

We  have  described  in  detail  a  thermal  fluid  electrochemistry 
model  for  the  study  of  the  electrochemical  performance  of  SOFC. 
The  model  couples  the  intricate  interdependency  among  the  ionic 
conduction,  electronic  conduction,  gas  transport,  electrochemical 
reaction,  thermal  conduction  and  heat  exchange  with  the  envi¬ 
ronment,  and  takes  into  account  the  contact  resistance  between  the 
cell  components  and  the  dependence  of  the  effective  electrode 
properties  on  the  microstructure  parameters  of  the  porous  elec¬ 
trodes.  The  validity  of  the  mathematical  model  is  demonstrated  by 
the  excellent  agreement  between  the  theoretical  and  experimental 
I-V  curves  of  mtSOFC  operated  at  different  temperatures.  The  nu¬ 
merical  model  is  used  to  study  the  effects  of  the  particle  sizes  and 
volume  fractions  of  Ni  and  LSM  as  well  as  the  anode  current  col¬ 
lection  mode  on  the  electrode  properties  and  current  production  of 
mtSOFC  with  the  typical  Ni-YSZ  anode,  YSZ  electrolyte  and  LSM- 
YSZ  cathode.  The  critical  importance  of  anode  current  collection  is 
revealed  and  the  current  output  of  mtSOFC  may  be  drastically 
increased  by  collecting  current  from  both  sides  instead  of  one  side 
of  the  anode.  This  finding  should  be  helpful  for  guiding  the  ex¬ 
perimentalists  to  design  improved  mtSOFCs. 

We  have  also  described  a  thermal  mechanical  model  for  the 
analysis  of  the  mechanical  behavior  of  SOFC,  with  the  thermal 
mechanical  properties  determined  by  the  available  theory  and 
experiment  and  linked  to  the  compositions  of  the  materials.  Based 
on  given  temperature  profiles,  e.g.,  the  results  by  the  simulations  of 
the  thermal  fluid  electrochemistry  model,  the  mechanical  model 
can  calculate  the  thermal  stress  distributions  of  the  mtSOFC  com¬ 
ponents,  followed  with  the  mechanical  failure  analysis  using  the 
Weibull  model.  Numerical  results  show  that  the  stresses  are  tensile 
for  the  Ni-YSZ  anode  and  mostly  compressive  for  the  YSZ  elec¬ 
trolyte  and  LSM— YSZ  cathode.  The  anode,  electrolyte  and  cathode 
suffer  higher  stresses  at  room  temperature  than  at  the  cell  oper¬ 
ating  temperature.  Flowever,  the  mechanical  failure  probability  is 
higher  at  the  cell  operating  temperature  than  at  room  temperature 
as  the  mechanical  strengths  of  the  materials  decrease  with  the 
increased  temperature. 

Optimizations  of  the  electrode  microstructure  and  composition 
for  both  the  electrochemical  and  mechanical  performances  of 
mtSOFC  are  considered.  Increasing  the  Ni  content  or  reducing  the 
Ni  particle  size  is  generally  helpful  for  improving  the  electro¬ 
chemical  performance.  Flowever,  the  tensile  stress  and  the 
accompanying  failure  probability  of  the  anode  and  the  compressive 
stresses  of  the  electrolyte  and  cathode  all  increase  with  the 
increased  Ni  content.  Overall,  a  Ni  volume  fraction  of  about  0.4  and 
a  Ni  particle  size  smaller  than  the  YSZ  particle  size  are 
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Nomenclature 

a},  aj ,  of  coefficient  of  the  specific  heat  of  material  species  i, 

(J  kg-1  I<-1) 

j bj,  bj,  bf  coefficient  of  the  thermal  conductivity  of  material 
species  i,  (W  m-1  K-1) 

Bo  permeability  coefficient,  (m2) 

Q  molar  concentration  of  species  i,  (mol  m-3) 

Cp  effective  specific  heat  of  the  fluid,  (J  mol-1  K-1) 

cref  Sutherland’s  constant. 

D  elasticity  matrix, 

Dz  diffusion  coefficient  of  species  i,  (m2  s-1) 

Dy(D?ff)  (effective)  binary  diffusion  coefficient,  (m2  s-1) 

Dfff  effective  Knudsen  diffusion  coefficient  of  species  i, 

(m2s-1) 

d  hydraulic  diameter,  (m) 

E  (Eeff)  (effective)  Young’s  modulus,  (GPa) 

Bo  Nernst  potential,  (V) 

Bh2(B02)  activation  energy  for  the  anode  (cathode),  (J  mol-1) 

G  (Geff)  (effective)  Shear  modulus,  (GPa) 

AG°  free  energy  change,  (J  mol-1) 

h  heat  transfer  coefficient,  (W  m-1  K-1) 

_  charge  transfer  rate,  (A  m-3) 
i  ei( 1  io)  vector  of  electronic  (ionic)  current  density,  (A  m-2) 
jo  exchange  transfer  current  per  unit  three  phase  boundary 

(TPB)  length,  (A  m-1) 

jg  ref (jg  ref)  anode  (cathode)  exchange  transfer  current  density  at 
the  reference  temperature,  (A  m-1) 

Jtpb  specific  transfer  current  density,  (A  m-1) 

I<i  bulk  modulus,  (GPa) 

ke ff  effective  thermal  conductivity,  (W  m-1  I<-1) 

kf  (ks)  effective  thermal  conductivity  of  the  fluid  (solid),  (W 
m-1  K-1) 

M  molecular  mass,  (kg  mol-1) 

Ni  molar  flux  of  species  i,  (mol  m-2  s-1) 

Nu  Nusselt  number 

n%  number  of  k-particles  per  unit  volume  in  a  composite 

electrode 

Pe d  (Pei)  Percolation  probability  of  electrode  (electrolyte)  particles 
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Pf  failure  probability  of  the  cell  component 

Ph2(Ph2o,  Po2)  partial  pressure  of  H2  (H20,  02),  (Pa) 

Q.  heat  source,  (W  m-3) 

Qact  (QentnQentr)  activation  (entropy,  ohmic)  heat  source,  (W  m-3) 
R  universal  gas  constant,  (J  mol-1  K-1) 

Ri  reaction  rate  of  species  i  at  TPB,  (mol-1  m-3  s-1) 

Pasr  contact  resistance,  (Q  cm2) 

rg  pore  radius,  (m) 

r/<  radius  of  k  particle,  (m) 

T  temperature,  (K) 

T0  room  temperature 

Text  external  temperature  of  air,  (I<) 

Tf  stress  free  temperature,  (K) 

TfW  temperature  of  the  furnace  wall,  (K) 

Top  cell  operating  temperature,  (K) 

7Vef  reference  temperature,  (I<) 

u  vector  of  convection  velocity,  (ms-1) 

u0  fuel  flow  rate,  (ms-1) 

Vceii  operating  cell  voltage,  (V) 

Z  average  coordination  number  of  all  solid  particles 

Zitk  average  number  of  /-particles  in  contact  with  a  fc-particle 


Greek  letters 

a  coefficient  of  thermal  expansion,  (K-1) 

a  com  coefficient  of  thermal  expansion  of  composite  material, 

(K-1) 


stress  tensor,  (MPa) 

initial  stress  or  the  residual  stress,  (MPa) 


el 

< 

^ten 


(of*)  effective  electronic  (ionic)  conductivity,  (S  m  ^ 


intrinsic  conductivity  of  material  k,  (S  m-1) 
tensile  stress,  (MPa) 

£rad  surface  emissivity 

?  strain 

?ei  elastic  strain 

?th  thermal  strain 

fo  initial  strain 

<pe i  {cpio)  local  electronic  (ionic)  potential,  (V) 
i pk  (V'f)  volume  fraction  of  /<- phase  (/-phase) 

\p[  percolated  volume  fraction  threshold  of  the  k  material 
rjact  activation  polarization,  (V) 

^act(^act)  anode  (cathode)  activation  polarization,  (V) 

^ohm^ohm)  electronic  ohmic  overpotential  of  anode  (cathode),  (V) 
a,i  /  c,i  eie,i  x  ionic  ohmic  overpotential  of  anode  (cathode, 
W  /ohm’  ohm;  electrolyte),  (V) 
tiasr  contact  ohmic  overpotential,  (V) 

VconiVcon)  concentration  balance  potential  of  anode  (cathode),  (V) 
p  density,  (kg  m-3) 

6  angle  of  particle  contact 

effective  volume  specific  TPB  length,  (m-2) 

effective  viscosity  coefficient  of  fluid  mixture,  (Pa  s) 

viscosity  coefficient  of  species  i  (j),  (Pa  s) 
viscosity  at  the  reference  temperature,  (Pa  s) 
tortuosity  factor 
(effective)  Poisson’s  ratio 


^TPB.eff 

M 

Mi  (Mj) 

Mref 

T 

u(ueff) 

V 


diffusion  volume,  (m3  mol  1) 
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